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Abstract 



Recent results from the SNO, Super-Kamiokande and Borexino experiments do 
not show the expected upturn of the energy spectrum of events (the ratio R = 
N bs/NssM) at low energies. At the same time, cosmological observations testify 
for possible existence of additional relativistic degrees of freedom in the early Uni- 
verse: AN e fj = 1 — 2. These facts strengthen the case of very light sterile neutrino, 
u s , with Amgi ~ (0.7 — 2) • 10 -5 eV 2 , which mixes weakly with the active neutrinos. 
The v s mixing in the mass eigenstate v\ characterized by sin 2 2a ~ 10 -3 can explain 
an absence of the upturn. The mixing of v s in the eigenstate 1/3 with sin 2 j3 ~ 0.1 
leads to production of v s via oscillations in the Universe and to additional contribu- 
tion AN e ff w 0.7 — 1 before the big bang nucleosynthesis and later. Such a mixing 
can be tested in forthcoming experiments with the atmospheric neutrinos as well as in 
future accelerator long baseline experiments. It has substantial impact on conversion 
of the supernova neutrinos. 



1 Introduction 



The large mixing angle (LMA) MSW solution [TJ [2] has been established as the solution 
of the solar neutrino problem [3], [I], [5], [6], [7], [8], [9], [10]. In assumption of the CPT 
conservation KamLAND confirms this result [TT], [12]. One of the main goals of further 
precision measurements of the solar neutrino fluxes is to search for possible deviations from 
the LMA predictions which would indicate new physics beyond the Standard Model with 
three mixed neutrinos. In particular, new physics can show up at the neutrino energies 
E = (1 — 7) MeV, i.e. in the transition region between the matter dominated conversion 
and vacuum oscillations. Here direct measurements of the spectrum are absent or inprecise 
and possible deviations from the LMA predictions can be relatively large. 

Some time ago in attempt to explain the low (about 2a) rate in the Homestake experi- 
ment [3j in comparison to the LMA expectation as well as the absence of clear low energy 
upturn of the spectra of events at SuperKamiokande and SNO we have proposed a scenario 
with light sterile neutrino, u s , which mixes weakly with active neutrinos [13]. Conversion 
of v e to v s driven by the mass squared difference Am^ ~ (0.2 — 2) ■ 10~ 5 eV 2 and mixing 
in the mass state V\ y sin 2 2a ~ 10~ 3 , leads to appearance of a dip in the v e — v e survival 
probability in the range (0.5 - 7) MeV which explains the data. 

After publication [13] several new experimental results have appeared which further 
support our proposal: 

• Measurements of the solar neutrino spectrum by SuperKamiokande-III [H] with lower 
threshold still do not show the upturn. 

• The SNO LETA analysis [13] gives even turn down of the spectrum in the two lowest 
energy bins. 

• The Borexino measurements of the boron neutrino spectrum also hint some tendency 
of the spectral turn down [16J. 

Although separately these results are not statistically significant, being combined they 
can be considered an evidence of some new sub- leading effect. 

At the same time, the cosmological observations indicate possible presence of additional 
radiation in the Universe in the epoch of last photon scattering. This is quantified by 
the effective number of neutrino species, N e ff, which is bigger than 3. Combined analy- 
sis of WMAP-7, measurements of BAO (Baryon Acoustic Oscillations) and new value of 
the Hubble constant H ) gives N eff = 4.34+g;|| [17]. WMAP-7 and Atacama Cosmology 
Telescope data lead to N e ff = 5.3 ± 1.3 (68 % C.L.) [18J. In the independent analy- 
sis [19] of these data the number of very light sterile neutrinos AN e ff = (0.02 — 2.2) (68 % 
C.L.) has been obtained. All this confirms the earlier finding based on the WMAP-3 data: 
Neff = 5.3±°;|±f;i±|;| [20]. 
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These results do not contradict the recent Big Bang Nucleosynthesis (BBN) bounds 
N e ff = 3.68^70 [21] (see discussion in [22] and theoretical considerations in [23]). Hence 
an additional radiation can be produced before the BBN epoch. 

In this connection we revisit our proposal of very light sterile neutrinos. We show that 
mixing of this neutrino in mass states v\ or/and z/2 can consistently improve description 
of the solar spectral data. We introduce mixing of this neutrino in the mass eigenstate v% 
which allows v s to be produced in the Early Universe with nearly equilibrium concentration, 
so that AN eff « 1. 

The paper is organized as follows. In sect. 2 we consider properties of the v e conversion 
in the presence of v s — mixing in the Sun generalizing our analysis in |13j . New feature, 
wiggles" in the survival probability, is described which appear for relatively large Amf at 
the E > 5 MeV. In sect. 3 we obtain bounds on the v s parameters from the Borexino 
measurements of the Be— neutrino flux. Spectra of the solar neutrino events have been 
computed for different experiments and confronted with the data. In sect. 4. the mixing 
of v s in V3 is introduced and phenomenological consequences of this mixing are studied, in 
particular, generation of v s in the Early Universe. The conclusion is given in sect. 5. In 
appendix we give some details of appearance of the wiggles in the survival probability. 



2 Sterile neutrino and conversion probabilities 

2.1 Generalities 

Let us consider the system of 4 neutrinos Vf = (u s , u e , z/ M , v T ) mixed in the mass eigenstates 
z/j, % = 0, 1,2,3. The sterile neutrino, v s , is mainly present in the mass eigenstate Vq with 
mass mo- It mixes weakly with active neutrinos and this mixing can be treated as small 
perturbation of the standard LMA structure. 

Coherence of all mass eigenstates is lost on the way to the Earth. Therefore the v e - 
survival probability at the surface of the Earth can be written as 

ft. = E 1-^21 W, (!) 

i 

where A s ei is the amplitude of the v e — > V{ transition inside the Sun and U ei = (v e \vi) is the 
element of the mixing matrix in vacuum. The quantities in eq. (pQ) satisfy the normalization 
conditions: J2i\U e i\ 2 = 1 and 

E \ A i\ = !• 

i 

During nights the solar neutrinos oscillate in the matter of the Earth. In this case U e i in 
eq. ([T]) should be substituted by the ^ — > v e oscillation probabilities inside the Earth, 
U ie ->• Af e , so that 

^ = EI4M 2 . (2) 
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In the production point the electron neutrino state can be represented in terms of the 
eigenstates in matter, u im , as 



$»im, (i = 0,l,2,3), 



where U™ is the mixing matrix element in matter in the production region. We denote by Aj 
the eigenvalues which correspond to the eigenstates v im . Introducing - the amplitudes 



of — > Vj transitions inside the Sun we can write 



X, E'7"-V (3) 

j 

Insertion this expression into (jSJ) gives 

^« = EIE^M a - ( 4 ) 



» j 



In the adiabatic case Ay = 5ij, so that 



^srapw- (5) 



For low energies we are interested in, the Earth matter effect is small and can be ne- 
glected in the first approximation. 

In what follows we will introduce mixing of the sterile neutrino in different mass eigen- 
states. In computations of effects for solar neutrinos we neglect the 1-3 mixing and there- 
fore consider the mixing of only three flavor states Uf = (i/ s ,i/ e , v a ) [y a is the mixture 
of and v T \ The mixing matrix which connects these states with the mass eigenstates, 
Vmass = (v , v \i v %)-, can be parametrized as 

U {3) = U e U a , (6) 

where U a is the matrix which mixes v s in v\ or/and V2 and Ug = 1/12(612) is the standard 
LMA mixing (the rotation by the angle 612 in the V\ — V2 plane). The Hamiltonian of the 
system in the v\ basis can be written as 

H f = U j e UlH dia9 U a U e + V, (7) 

where 

H diaa = diag(^ , H u H 2 ) = ^diag(m2, m?, ml), 

V = diag(0,K,K) (8) 

are the diagonal matrices of the eigenvalues of the Hamiltonian in vacuum and the matter 
potentials correspondingly; V e = V^Gjr(n e — 0.5n n ) and V a = —(l/y/2) Gpn n , with n e and 
n n being the electron and neutron number densities. 
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Figure 1: The level crossing scheme for the neutrino energy E = 10 MeV. Dependence of the 
eigenvalues of total Hamiltonian in matter, Aj, on distance from the center of the Sun. The 
LMA neutrino oscillation parameters are taken as Am^ = 8 x 10~ 5 eV 2 and tan 2 # 12 = 0.44, 
while sterile neutrino parameters are i?A = 0.25 and sin 2 2a = 10~ 3 . 



It is convenient to consider effects of sterile neutrino mixing in the basis rotated by the 
1-2 mixing in matter U$ m \ (u s , v^ A , vkrn A \ which would diagonalize the Hamiltonian in the 
absence of mixing with sterile neutrinos (i.e. when U a — I). In this basis the Hamiltonian 
becomes 

H a = U 6m H f Ui m = U em U}ulH d ^U a U e Ul + U 6m Vll\ m . (9) 
Since U a is small rotation we can represent it as 

U a = I + U s , (10) 

where Us = U a — I. Inserting this expression into j9]) and taking the lowest order terms in 
Us we obtain 

H a = H^ + H s , (11) 

where 

H*" = diag(tf , Af MA , \ L 2 MA ) (12) 
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is the Hamiltonian in the absence of mixing with sterile neutrino and 

H s = Ui e U}H dia °U Ae + h.c. 



(13) 



is the correction to the Hamiltonian due to mixing with sterile neutrino. Here Ad = 
and 

/ 1 



Uab = UoU, 



t 



cos(fl-# m ) 
V sin(# - e m ) 



— sin(# — Q Tl 
cos(8 — dm 



(14) 



We denote the ratio of mass squared differences as 

Ra = Am oi 

Depending on mass and mixing of the sterile neutrino (i.e. the form of U a ) one can obtain 
several phenomenologically different possibilities. 



2.2 The case m\ < tuq < mi 

This case corresponds to *C 1 (as in [13J). For the neutrino with energy E = 10 MeV 
the level crossing scheme which gives dependence of Aj, (i = 0, 1, 2) on the distance inside 
the Sun (or on the density), is shown in fig. [TJ Details of construction of this scheme can 
be found in [13]. For other energies the scheme can be obtained from the one in Fig. CD by 
shifting the picture with respect to the frame to the right with increase of energy and to the 
left with decrease of energy. According to fig. [T]for E = 10 MeV the sterile neutrino level, 
\ s , has two resonances - two crossings with the original (without u s ) level Af M . at smaller 
density rzf , and at higher density, n^, {n^ > nf). With increase of Am^ the density 
increases, whereas - decreases, they approach each other and then merge. 

The flavor content of the eigenstates in matter, z/ im , i.e., the mixing matrix elements 
f/™(n(r)) as functions of density (distance from the center of the Sun) is shown in fig |2J 
Notice that since the v s mixing in the is absent, the change of flavor of V2 m with density 
is the same as in the LMA case: = U^ LMA . The intersections of the lines correspond 
to resonances. 

Consider evolution of this system from the production point in the central region of the 
Sun to the solar surface. The eigenstates z/ 2m and z/ 3m evolve adiabatically 

i.e., A 2 j = A 3 j = 5%j and therefore according to ([3]) 

A S e2 = U™, A s e3 = US ^ U e3 . (15) 
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Figure 2: Flavor content of the mass eigenstates in matter for the same neutrino parameters 
as in fig. [TJ 



In the last equality we have taken into account that due to large value of Am^ the 1-3 mixing 
is practically unaffected by the solar matter in the energy range of interest. Inserting the 
amplitudes (|T5j) into (JTJ) we obtain 

P ee = L4fi| 2 |[/ el | 2 + |Af | 2 |f/ e0 | 2 + \U2\ 2 \U £ 2\ 2 + \U e3 \ 4 . (16) 

In turn, the amplitudes of v e — transitions to v\ and to uq can be written according to ([31) as 

A s el = U™A n + U™A 01 , Af = U™A 10 + U™A 00 . (17) 

Finally, insertion of ( TTTI) into ( fTBT) gives 

P ee = \U™A n + U™A m \ 2 \U el \ 2 + \U2A 10 + U™A 00 \ 2 \U e0 \ 2 + \U™\ 2 \U e2 \ 2 + \U e3 \\ (18) 

Since \U e0 \ 2 ~ 10 -3 and \U e3 \ 4 ~ 4 ■ 10~ 4 , the probability ([TBI can be written approximately 

as 

P ee w \U3A n + U™A m \ 2 \U el \ 2 + m\ 2 \U e2 \ 2 . (19) 

With decrease of neutrino energy effectively the pattern in fig. [T] shifts to the left. 
Therefore at low energies there is no sterile neutrino resonances, the evolution proceeds 



adiabatically and 

P ee ^\U2\ 2 \U el \ 2 + \U2\ 2 \U e2 \ 2 . (20) 

Furthermore, E/™ ~ JJ^ LMA and U'£, ~ U^ LMA , where U™ LMA are mixing parameters in 
the pure LMA (2z/) case without sterile neutrinos. So, P ee ~ P^ e MA . 

With increase of energy the low density resonance becomes effective (the corresponding 
level crossing scheme for E = 8 MeV is shown in fig. 1 of our paper p2]). In the adiabatic 
case (relatively large U i) we have A 01 w 0, An ~ 1 and the expression for P ee in eq. f fT9|) 
is reduced to the one in (120]) . However now [7™ ~ if neutrino is produced above the 
resonance layer, and consequently, P ee ~ \U^ \ 2 \U e 2\ 2 . 

In the non-adiabatic case: 

p ee « i^ 01 | 2 if/ el | 2 + ra 2 it4 2 i 2 ~ i^ LMA A 01 | 2 |f/ el | 2 + i^i 2 if/ e2 | 2 , (21) 

where we have taken into account that U'^ ~ U™ LMA . In the case of strong adiabaticity 
violation (for very small Uqi) when Aoi ~ 1, the probability in ( 12T|) is reduced to the 
standard LMA probability. 

With further increase of energy at E > 9 MeV two sterile resonances are realized and the 
amplitudes Aqi and An can be written in terms of transition amplitudes in each resonance 
A$ (a = 1,2) as 

a - A {2) A {1) + A {2) A {1) - A {2) A {1) + A (2) A (1) (22) 

Sill — ^11 ^11 T ^10 ^01 ) ^01 — ^00 ^01 + ^01 ^11 • 

One can get different results depending on the adiabaticity conditions in each resonance. If 
the crossings are adiabatic, A i w 0, An « 1, we obtain from ( 1191) 

m\2\TT 1 2 i I r rm 1 2 I 



p ee « \u2nu e i\'+\u2\ z \u e2 \'. 

This expression is similar to the one for the standard LMA case, however now \U™\ 2 = 
1 - \U™\ 2 - \U™\ 2 = \U™ LMA \ 2 - TO 2 , and consequently, 

P — pLMA _\tt \2\jjm\2 
r ee — *ee \ U el\ \ u e0 I ■ 

Here the second term describes the dip in the adiabatic case. (Notice that due to smallness 
of v s mixing in vacuum U e \ ~ U^{ MA .) 

If adiabaticity is strongly broken in both resonances due to smallness of sterile mixing, 
then A$ « A$ « 1 and A$ « a£ 3 « « 0. Therefore according to (|22jl, Am » 0, 
An ~ 1, and as in the adiabatic case: 

P ee ^ |^| 2 |f/ el | 2 + |^| 2 |f/ e2 | 2 . (23) 

However, now the mixing parameters are approximately equal to the standard LMA param- 
eters without sterile neutrinos. Therefore P ee ~ P^ e MA . 
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In the case of mixing of v s in v\ and vq only, the matrix U a equals 

/ cos a sin a \ 
Z/qj = I — sin a cos a . (24) 




Explicitly the flavor mixing can be parameterized as 

Vq = cos a v s + sina(cos6 l 1 2 v e — sin 9 12 v a ), 
ui = cos a (cos6 l 1 2 v e — sin 6^2 v a ) — sin a u s , 
v 2 = sin 9 12 v e + cos 9\ 2 v a . 

Here u e and u a (a combination of and u T ) mix with the angle 9± 2 in the mass eigenstates 
V\ and v 2 having the mass split Amf 2 - In terms of these mixing angles: 

U e \ = cos a cos #12, f/ e o = sin a cos 6*12, U e2 = sm9\ 2 . 

With this parametrization eqs. ( fl9|) and (|23|) reproduce the corresponding results of our 
previous paper [13] . 

From eqs. (TTT1 - we find using (124]) that the Hamiltonian in the rotated basis equals 

(0 — A ^ 01 sin o; cos(6 > — 9 m ) smasm(9 — 9 m ) \ 

\f MA -H ' • (25) 

\ L 2 MA -H Q ) 

Due to smallness of a the off-diagonal terms are much smaller than the diagonal ones. If 
Ra C 1, so that H is close to \[ MA and crosses this level, there is no crossing of \ 2 MA and 
H Q levels, and the state v 2 ^ A decouples. Then mixing of v s and v^ A is determined by 

sin a— ^cos(# - 9 m ) = sma—^-R A cos{9 - 9 m ). (26) 

Since the transition occurs in the resonance region, 9 m should be taken at the density which 
corresponds to the sterile neutrino resonance. The expression for mixing (|26|) allows to 
understand behavior of the conversion probability on Ra, a and neutrino energy. 

In figs. [31 H]we show results of numerical computations of the v e survival probability P ee 
as functions of the neutrino energy. We have performed complete integration of the evolution 
equations for the 3^-system and also made averaging over the production region of the Sun. 
The following analytical consideration allows to understand the numerical results shown in 
fig. [3j The v s — mixing and v e — > v s conversion lead to appearance of a dip in the energy 
dependence of the v e — survival probability. The survival probability in the pure LMA case 
is given by 

Pee MA ~ \U™ LMA \ 2 \U^ MA \ 2 + \U™ LMA \ 2 \U^ 2 MA \ 2 + |f/ e3 | 4 . (27) 
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Figure 3: The survival probability of the electron neutrino as function of neutrino energy 
for different values of the sterile-active mixing parameter sin 2 2a and mass scale <C 1 . 
The active neutrino parameters are Am^ = 8 x eV 2 and tan 2 9 = 0.44. 




Since u s mixing in the v 2 is absent, U^ 2 MA = U e2 and U^ LMA = and the probability ( 1271) 
can be rewritten as 

Pee MA ~ l^eTH^ell 2 + + \U e3 \ 4 . (28) 

Using normalization conditions 

i^ii 2 + i^oi 2 + ra 2 = i 

(we neglect 1-3 mixing here) and J2i \U e i\ 2 = 1 (i = 0,1,2,3) as well as the expressions in eqs. 
( fl6l) and ( 1281) we can find difference of the probabilities with and without sterile neutrino 
effect: 

AP ee = Pte MA -Pee = \Al Q \\l-\UM-\Ue0\\2\A^-l-\U-\ 2 ) 

w |Af | 2 (l-|[/ e2 | 2 )^P es (l-|[/ e2 | 2 ), 

where P es ~ |^4fol 2 i s the probability of v e — > u s transition. The quantity AP ee describes 
the dip which has the following properties (see fig. [3] and also discussion in [T3]): 
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Figure 4: The same as in fig. [3] for higher values of mixing angle a. 



1. A position of the dip (its low energy edge) is given by the low density resonance taken 
at the central density of the Sun Ei(n c ). With increase of Am^ the dip shifts to higher 
energies. 

2. Maximal suppression in the dip depends on and a. For small i?A (large spit 
between the two resonances) and large a (sin 2 2a > 1CT 3 ) the absolute minimum can be 
achieved at the adiabatic crossing of the v s — resonances. With increase of (smaller split 
of the resonances) or/and decrease of a (stronger violation of the adiabaticity) a suppression 
in the dip weakens. Also with decrease of a the dip becomes narrower. 

3. For large Am^ and relatively small a the survival probability as function of the 
neutrino energy has wiggles (see fig. [3]). The wiggles are result of interference of the 
two amplitudes in the first term of ( fT9l) which develops over finite space interval. Indeed, 
according to ( fl~9]) there are two channels of transition of v e to V\ : 

(i) v e has admixture U™ in v lm , the latter adiabatically evolves to V\\ v e — > V\ m —> u±, 
and the amplitude equals U™An. 

(ii) u e has admixture in uo m ; this state transforms to V\ due to non-adiabatic tran- 
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sition: v e — > z/ 0m — > v\. The corresponding amplitude is U^A i. 

The two contributions to the amplitude interfere leading to the oscillatory dependence 
of the probability on energy (wiggles). Introducing P 01 = |A i| 2 , so that \A U \ 2 = 1 — P i> 
we can rewrite the probability (|T9|) as 



P ~ 



u el \ 



U™\ 2 (1 - Poi) + m\ 2 P 01 + U2U% cos0 J P 01 (l - P, 



01 j 



\U, 



< 2 



\2\jjm\2 
\ U e2\ i 



where <p = arg(AQ 1 A n ) and we assumed for simplicity that U™ and t/^J are real. The 
oscillatory behavior follows from the energy dependence of the phase <fi. The key point is 
that the phase is collected over restricted space interval, L, and therefore is not averaged 
out even after integration over the production region. Indeed, the phase is acquired from 
the neutrino production point to the second (low density) resonance. Below the second 
resonance (in density) both "trajectories" (channels of transition) coincide. Appearance of 
the wiggles requires the adiabaticity violation. In the adiabatic case A m = and only one 
channel exists. Unfortunately, it will be difficult, if possible, to observe experimentally these 
wiggles. Some more details concerning the wiggles are presented in the Appendix. 



If v s mixes in i>2, then 

/ cos a' sin a' 

U a = 10 

y — sin a' cos a' 

and the Hamiltonian can be obtained from fT23|) by substitutions 

cos(# — 9 m ) — > sin(6 l — 9 m ), sin(6 l — 9 m ) —> — cos(6 ) — 9 m ), Am^ — > Am 



02- 



(29) 



(30) 



Again the state v 2m decouples and v s — mixing is given by 



sin a 



sm(6> - 9 m ) = sin a - 



2P 



2P 



1 - P A ) sin(# - 9 n 



(31) 



Notice that this mixing appears due to matter effect and it is absent in vacuum when 9 m — > 9. 
It happens that for values of Pa we are considering (1 — Pa) sin(# — 9 rn ) ss Pa cos(9 — 9 m ) 
and therefore the probabilities in this case are very similar to those shown in figs. [31 HI 



2.3 mo > vrt2 > mi and other possibilities 

For mo > rri2 > rn\ the ratio Pa > 1. Since below the LMA resonance A2 = \\ MA and 
V n have practically the same dependences on density (radius) (see fig. [TJ), there is only one 
crossing of X s with \z MA , and there is no crossings for Ai?1q 2 < 0. Now the evolution of 
states v lm and z/ 3m is adiabatic, so that 

aS ^ jjm _ TjOm aS ^ tj 
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Consequently, 



P ee = \U2\ 2 \U el \ 2 + \A s e2 \ 2 \U e2 \ 2 + \A s J 2 \U eQ \ z + |[/ e3 | 
where 

A s e2 = U™A 22 + C/-A 02 , Af Q = E^Aa, + US^oo- 



These expressions are similar to the expressions in (1161) and ([3]) with interchange of indexes 
1 o 2. 

In the adiabatic case we have 

P ee = \U2\ 2 \U e i\ 2 + \U%\ 2 \U e2 \ 2 + \U™\ 2 \U e0 \ 2 + |[/ e3 | 4 . 

Now the effect of sterile neutrino is due to difference of and U^ LMA : \U^\ 2 = \U^ LMA \ 2 - 
\U'^\ 2 . In the strongly non-adiabatic case A 20 « A 2 ~ 1 and 

P ee = \U2\ 2 \U el \ 2 + \U™\ 2 \U e0 \ 2 + \U™\ 2 \U e2 \ 2 + \U e3 \\ (32) 

with [7™ « 0, and C/^ « U^ LMA , so that ([32]) is reduced to the LMA probability. 

For v s — mixing in Vq and V\ the Hamiltonian is given by the same expression as in eq. 
( 1251) . However, now H crosses X 2 MA and the state v\^ A decouples. According to ( |25j) . the 
mixing of v s and is determined by 

sin a 01 sin(6> - 9 m ) = sin a 21 R A sm(9 - 9 m ), (33) 

and this mixing is due to matter effect. For v s — mixing in vq and v 2 performing substitutions 
f l30l) we obtain the z/ s — v 2 ^ A mixing 

sin a'— -52 cos(0 - m ) = sin a'— -21(1 - i? A ) cos(# - 9 m ). (34) 
2_c/ 2i!/ 

In the first case, eq. (1331) . the mixing oc i?A sin(6> — 9 m ) is larger than in second case ( 1341) : 
oc (1 — -Ra) cos(# — # m ), since i?A ~ 1- Furthermore, the first mixing increases with energy: 
the sterile resonance is above the LMA resonance and therefore 9 m > 45°; this angle, and 
consequently | sm(9 — 9 m )\, increase. As a result, the effect does not disappear at high 
energies (see and ED • 

In comparison to the case i?A "C 1, now v e has smaller admixture in v 2 , \U e2 \ < \U e i\, 
however the initial admixture of v e in v 2m can be larger: |C/^| > Therefore the 

overall effect is large (see figs [5] and [6]) . Indeed here we have only one level crossing and 
improvement of the adiabaticity in the resonance leads to stronger transition. With increase 
of Amg 2 and therefore, the dip moves to high energies but the resonance shifts to higher 
densities, i.e., to the central regions of the Sun where the density gradient is smaller and 
adiabaticity is better. Here substantial change of the probability exists for smaller mixing 
angles. 
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Figure 5: The same as in fig. |3]for _Ra > 1- 



If m < mi < m 2, so that Am^ < 0, the sterile level X s crosses v\^ A at high densities 



only 



The resonance energy equals 



n^2n LMA (l + R A ). 



E 



R n ' 



where n n is the number density of neutrons. In this case we have the same general ex- 
pressions for the survival probability as in ([15]) and f|T9|) . Consequently, the expressions for 
P ee in adiabatic and non-adiabatic limits coincide with those in e.g. fj2ip for one sterile 
resonance. However the dip here is at high energies. 



In the case of flavor mixing, that is the mixing of u s with v ei v a the matrices Uq and U a 
should be permuted, so that = U a Ue (compare with ([6])). It can be shown that now 
the off-diagonal elements of the Hamiltonian H a contain terms with Am^ and Am^ simul- 
taneously. As a result, the probabilities have energy dependences which are intermediate 
between those for mixings in mass states v\ and v 2 . 
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Figure 6: The same as in fig. H]for i?A > 1 and v s mixing in v^. 

3 Solar neutrino data and sterile neutrino effect 



In what follows we will consider scenario with m\ < mo < mi. This possibility gives better 
description of the data: it leads to significant modification of the survival probability in the 
transition region and weakly affects spectra at high energies. 



3.1 Borexino measurements of the Be-neutrino line 

The results of Borexino experiment [10] are in very good agreement with prediction based 
on the LMA solution and the Standard Solar Model. Within the error bars no additional 
suppression of the flux has been found on the top of P^ e MA . In Borexino (and other ex- 
periments based on the ^e-scattering) the ratio of the numbers of events with and without 
conversion can be written as 

RBorexino = PeeiX Tj-\-T vP es , (35) 
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where r = a{y^e — v^e) / a{y e e — u e e) is the ratio of cross-sections. Using Eq. (1351) we find 
an additional suppression of the Borexino rate in comparison with the pure LMA case [13] : 

^R-Borexino = Rsorexino ~ ^Borexino = (1 ~~ r )^Pee + r Pes ~ A.P ee (l r tan #12). 

In fig. Owe show dependence of the survival probability at E = Ese as function of for two 
different values of the mixing angle a. We show the Borexino bounds on P ee obtained from 
the experimental result [IE] and relation (1351) . According to this figure for sin 2 2a = 10 -3 
the range Ra = 0.005 — 0.072 is excluded at la level. For sin 2 2a = 5 ■ 10~ 3 we obtain 
slightly larger exclusion interval: i?A = 0.001 — 0.075. The Be— neutrino line can not be in 
the dip or the dip should be shallow which then will have little impact on the higher energy 
spectrum. So, essentially the allowed values of masses (which influence the upturn) are 

R A > 0.075, or Am 2 01 > 0.5 • 10~ 5 eV 2 . 




Figure 7: The survival probabilities for the monoenergetic 7 Be— (red) and pep— (black) 
neutrino fluxes as functions of R&, for two values of mixing: sin 2 2a = 1 x 10 -3 (solid lines) 
and 5 x 10~ 3 (dashed lines). The active neutrino oscillation parameters are the same as in 
fig. [3] The horizontal line and shadowed band show the central value and la band for the 
suppression factor determined by Borexino. 
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3.2 Upturn of the boron neutrino spectrum 

Using the survival probabilities obtained in sect. 2 we have computed the energy spectra of 
events for different experiments with and without sterile neutrino. These spectra together 
with experimental data are presented in figs. M - CE21 We did not searched for the best 
fit of the data points, and the figures have an illustrative character. Notice that due to 
uncertainty in the original boron neutrino flux the experimental points can be shifted with 
respect to the theoretical lines by about 15%. 




10 12 

E - (MeV) 



Figure 8: Prediction for the SuperKamiokande-I spectrum versus experimental data [6]. 
The active neutrino parameters are the same as before; the sterile neutrino parameters 
equal = 0.20 and sin 2 2a = 1 x 10~ 3 (red line) and 5 x 10~ 3 (blue line). The pure LMA 
spectrum is presented by the dashed black line, with a normalization factor / = 0.91 to 
reproduce the total observed number of events. We use the 8 B— neutrino flux according to 
the GS98 solar model [21]. 



In figs. Eland [9] we show the ratio of the number of events in SuperKamiokande-I (SK-I) 
with and without oscillations for two different values of Ra- Different curves correspond to 
the standard LMA solution (dashed) and the spectra with conversion to sterile neutrino. In 
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Figure 9: The same as in fig. E]for = 0.08. 



the presence of sterile neutrino mixing the upturn can be completely eliminated and even 
transformed into turn down of the spectrum. In fig. [8] the dip at E ~ 4 MeV corresponds 
to the dip in the probability at approximately the same energy as in the fig. |3] (middle 
pannel). The difference of the predictions with and without sterile neutrino can be as big 
as (15 - 20)% at E e < 5 MeV. 

The SuperKamiokande-III data (SK-III) (fig. [T01 has additional lower energy [T4], how- 
ever statistics is lower than in SK-I. Again there is no clear indication of the upturn in the 
SK-III spectrum and theoretical lines with sterile neutrino mixing can describe the data 
better than pure LMA solution. 

SNO (fig. [TTl) is more sensitive to distortion of the neutrino spectrum. However, the 
dip in the electron spectrum is shifted to low energies by the threshold of the CC reaction 
on the deuteron: E = 1.44 MeV. Experimental points are from the SNO-LETA charge 
current event analysis [12] . Two low energy points of the spectrum show a sharp turn down. 
This can not be reproduced by the proposed dip, although with the dip the description is 
better 0. Also the Borexino spectrum (fig. [12]) can be fitted better in presence of sterile 
neutrino mixing. 

1 Too sharp decrease of signal in the lowest energy bins is probably statistical fluctuations or some 
systematics. 
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Figure 10: Prediction for SuperKamiokande-III spectrum versus experimental data [14J. 
The neutrino parameters and the solar model as well as the normalization factor for pure 
LMA spectrum are the same as in fig. [8] (left). 

According to fig. |E]-[D±]an improved description of the data can be achieved with 
Am 2 01 > 1.5 ■ l(T 5 eV 2 , sin 2 2a ~ 1(T 3 . 

3.3 Further tests 

With higher statistics Borexino will improve precision of measurements of the boron neutrino 
spectrum. Also SuperKamiokande will achieve better measurements of spectrum at lower 
energies. The KamLAND solar (23] and SNO+ experiments will further check the 
presence of the dip. 

Additional probe of the existence of sterile neutrino (and restriction on its parameters) 
can be provided by measurements of the pep— neutrino line with E = 1.44 MeV since 
the pep— neutrino flux is known with high precision. In fig. [7] we show dependence of the 
suppression factor for the pep— neutrinos as function of i?A- With increase of i?A the dip 
shifts to higher energies. In the interval = 0.07 — 0.11, the Be— neutrino flux has the 
LMA suppression, whereas the pep— flux can be suppressed by factor 0.15 - 0.20 (the LMA 
suppression is 0.52). In the range > 0.12 both fluxes have the LMA suppression. 
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Figure 11: Prediction for the SNO-LETA electron spectrum versus experimental data [15]. 
The neutrino parameters and solar model are the same as in fig. [HJ 



In the range Ra > 0.12 the CNO— neutrinos are sensitive to the presence of the dip, 
however the original fluxes of these neutrinos are not well known. 

4 Extra radiation in the Universe and v s — z/3 mixing 

Smallness of mixing of the sterile neutrino in the states v\ or/and ^2 (|£4i| 2 < 10 -3 ) does 
not lead to production of significant concentration of v s in the Early Universe via neutrino 
oscillations [13J. However, substantial abundance of v s can be generated if v s mixes in the 
state z/ 3 and is large enough. Description of the solar neutrino data presented in the 
previous sections does not change substantially, if v s mixes with combination 

v' T = cos 6*23 v T + sin6» 23 z/ M ~ ^3, 

where #23 is the standard 2-3 mixing angle. The v s — v' T mixing can be parametrized by 
the angle f3 as 

z/ 3 w cos fiv T + sin [3v S) z/ ~ cos fiv s — sin fiv' T , (36) 

so that U S 3 ~ sin/3. Here we neglect small rotations by the angles a and #13 which do not 
influence conclusions of this section. (These mixings can be introduced before or after the 
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Figure 12: Prediction for B— neutrino spectrum at Borexino versus with experimental data 
[16] . The neutrino parameters and solar model are the same as in fig. [HJ 



rotation fl36|) ). Since Am^ <C Am^ <C Am^, the mass squared difference of u 3 and z/ 
equals 

Am 2 « Am^ = 2.5 ■ eV 2 . 
For this value of Am| the mixing angle ft is restricted by the atmospheric neutrino data 

sin 2 p< 0.2 -0.3, (90% C.L.) 

and by the MINOS searches for depletion of the neutral current events [28]. For zero 1-3 
mixing the bound ft < 28.8° has been established [28] which corresponds to 

sin 2 ft< 0.23, (90% C.L). 

In the presence of non-zero 1-3 mixing the bound becomes much weaker. 

If sin 2 (3 ~ 0.2, then according to [29] the sterile neutrinos practically equilibrate before 
the BBN epoch both in the resonance channel and in non-resonance channels, i. e. in neutrino 
and antineutrino channels. Consequently, in the epoch of nucleosynthesis and latter the 
additional effective number of neutrinos is 

AN eff « 1. 
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The value AiV e // ~ 0.8 can be obtained for sin 2 /3 ~ 0.03 in the non-resonance channel and 
sin 2 (3 ~ 10~ 3 in the resonance channel. According to [27] AN e ff « 0.8 is generated, if the 
v s — mixing is about sin 2 (3 = 0.02. 

The CNGS experiment has also some potential to restrict sin 2 (3 [30\. 

Let us consider other phenomenological consequences of the v' T — v s mixing. The level 
crossing scheme can be obtained from fig. [1] by adding the third active neutrino level and 
expanding whole the picture to the left. With increase of density the A2 increase until the the 
region of 1-3 resonance and then turns down and decreases in parallel to Ai. Consequently, 
the sterile level X s ~ Ao (horizontal line) will cross A2 at some density above the 1-3 resonance 
density. Thus the mixing of v s in u 3 leads to appearance of the resonance in v' T — v s channel 
(normal mass hierarchy) at the density determined by 

1 „ Am 2 ,, Ami, 

V n = -^G F n n « — « — 
^2 IE 2E 

In the isotopically neutral medium this density is about 2 time larger than the density of 
1-3 resonance. For the inverted mass hierarchy the resonance appears in the antineutrino 
channel v' T — v s . 

Inside the Earth the v' T — v s resonance energy equals E « 12 GeV and wide resonance 
peak appears in the range (10 - 15) GeV. This can be tested in studies of the atmospheric 
neutrinos (spectra, zenith angle dependences) in the IceCube DeepCore detector [31] and 
in next generation Megaton-scale experiments [32] . Effect of such a mixing should show up 
in the long baseline experiments as the energy dependent disappearance of the v^— flux. 

The v' T — v s mixing also influences the supernova (SN) neutrino conversion. The corre- 
sponding level crossing in the collapsing star will be adiabatic (at least before shock wave 
arrival) and therefore v' T converts almost completely in this resonance into v s . At larger 
distances from center of a star this v s — flux will encounter the lower density v s resonances 
due to the crossing of v s and v\ m levels (see fig. H]). The latter will lead to partial conversion 
of v s into z/ e , since the adiabaticity is broken in these resonances. Hence, the following chain 
of transitions is realized: 

K( v i» v t) v s -> v„, v e . (37) 

Consequently, even for relatively large 1-3 mixing which leads to the transition u e — > z/3 with 
I (^el^) 1 2 *C 1 (for normal mass hierarchy), the v e signal may not be strongly suppressed due 
to conversion described in (137|) . In the case of inverted mass hierarchy similar consideration 
holds for the antineutrino channels. 

In this consideration for simplicity we have neglected possible collective effects due to 
neutrino-neutrino scattering and effects of shock wave propagation (see [32J). 
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5 Conclusions 



1. Recent measurements of the energy spectra of the solar neutrino events at SuperKamiokande, 
SNO, Borexino do not shown the expected (according to LMA) upturns at low energies. 
The absence of the upturn can be explained by mixing of very light sterile neutrino in the 
mass states v\ or/and z/ 2 with Am^ ~ (0.7 — 2) • 1CT 5 eV 2 (Ra = 0.07 — 0.25) and mixing 
sin 2 2a = (1 — 5) ■ 10~ 3 . Such a mixing leads to appearance of the dip in the v e — survival 
probability in the energy range (1-7) MeV, thus removing the upturn of the spectra. For 
Am§j ~ 2 ■ 10~ 5 eV 2 and sin 2 2a ~ 5 • 10~ 3 the u e — u s conversion can even produce a 
turn down of the spectra. Description of the existing solar neutrino data in the presence of 
mixing with sterile neutrino is apparently improved. 

Values of Am 2 < 0.6 ■ 10~ 5 eV 2 (for mixing angle interval sin 2 2a = (1 - 5) ■ 10 -3 ) 
are excluded by the Borexino measurements of the Be— neutrino flux. The presence of the 
dip can be further tested in future precision measurements of the low energy part of the 
B— neutrino spectrum as well as the pep— neutrino flux. 

2. Mixing of u s in the u 3 mass eigenstate with \U s s\ 2 ~ 0.02 — 0.2 leads to production of 
significant concentration of v s via oscillations in the Early Universe. For \U s z? ~ 0.1 — 0.2 
nearly equilibrium concentration can be obtained both in neutrino and antineutrino chan- 
nels thus generating additional effective number of neutrinos AN e jf ~ 1 before the BBN 
epoch. This can explain recent cosmological observations. 

3. Mixing of v s in z/3 leads to a number of phenomenological consequences, in particular, 
it can affect the atmospheric and accelerators neutrino fluxes as well as fluxes of the SN 
neutrinos. The mixing leads to existence of the v s — v' T resonance. For neutrinos crossing 
the Earth the resonance should appear at energies E ~ 10 — 15 GeV. This can be tested 
in future atmospheric neutrino studies with Megaton-scale detectors as well in the long 
baseline experiments with accelerator neutrino beams. 
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Appendix A: Wiggles 

As we described in sect. 2, wiggles in the dependence of the v e — survival probability on 
energy are the result of interference of the amplitudes which contribute to the same u e — > v\ 
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transition. The zoomed view of the survival, P ee , and transition, P es , probabilities is shown 
in fig. [13j The period of wiggles is about (0.5 - 0.6) MeV. 
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neutrino energy 



Figure 13: A zoomed view of the survival and transition probabilities in the energy range 
where the wiggles can be well seen. The sterile neutrino parameters equal = 0.25 and 
sin 2 2a = 10~ 3 . 

The key feature which leads to the wiggles with rather large period in the energy scale 
and therefore prevents them from being averaged out at the integration over the production 
region is that the interference phase is collected over relatively small distances L. These 
are the distances between the production point and the low density v s — resonance or the 
distance between the two v s resonances as can be seen in the fig. O For the neutrino energy 
E ~ 8 MeV the distance L w 20/ m where l m is the oscillation length in matter. Therefore 
the period of wiggles can be estimated as AE/E ~ l m /L ~ 1/20 in agreement with results 
of fig. M 



The wiggles are partially averaged due to integration over the production region. Notice 
that with decrease of Am^ the lower resonance shifts to lower densities and the distance 
L increases leading to smaller period of wiggles and stronger averaging. This is one of the 
reasons of disappearance of wiggles with decrease of Am^. The amplitude of wiggles also 
decreases with increase of a: the latter means better adiabaticity and therefore suppression 
of the contribution of one of the channels responsible for interference. 
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Figure 14: Transition probability to sterile neutrino for neutrino created at the Sun center 
as function of distance from the center of the Sun. The sterile neutrino parameters are 
R A = 0.25 and sin 2 2a = 10~ 3 . 
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